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ABSTRACT 


Refined  mathematical  methods  are  required  for  the  ana¬ 
lytical  solution  of  the  partial  differential  equation  governing 
steady,  two-dimensional,  compressible,  transonic,  potential 
fluid  flow.  This  equation  is  nonlinear  in  the  physical  plane 
and  so  does  not  lend  itself  to  standard  analytical  solution 
methods.  The  Molenbroek-Chaplygin  transformation,  where  the 
physical  Cartesian  coordinates  as  the  independent  variables 
are  replaced  by  the  velocity  magnitude  and  direction  as  the 
independent  variables,  linearizes  the  governing  equation  which 
may  then  be  analytically  solved.  The  plane  where  the  said 
velocity  parameters  are  the  independent  variables  is  termed 
the  hodograph  plane.  Likewise,  the  transformed  differential 
equation  is  known  as  the  hodograph  equation  and  it  is  solved 
by  hodograph  methods. 

This  mathematical  study  addresses  the  solution  of  tran¬ 
sonic  flow  phenomena  by  an  extension  of  Lighthill's  hodograph 
method.  Lighthill's  method  transforms  a  given  solution  of 
the  Laplace  equation  into  a  solution  of  the  hodograph  equation 
for  subsonic  flows  only.  A  new  relation  is  developed  in  this 
study  extending  this  transformation  technique  to  include  flows 
up  to  Mach  2.2735  in  air.  Requiring  only  numerical  data  con¬ 
cerning  the  velocity  field,  this  hodograph  method  is  compu¬ 
tationally  efficient  and  mathematically  straightforward. 
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ABSTRACT 


Refined  mathematical  methods  are  required  for  the  analytical 
solution  of  the  partial  differential  equation  governing  steady,  two- 
dimensional,  compresssible,  transonic,  potential  fluid  flow.  This 
equation  is  nonlinear  in  the  physical  plane  and  so  does  not  submit 
itself  to  standard  analytical  solution  methods.  The  Molenbroek- 
Chaplygin  transformation,  where  the  physical  Cartesian  coordinates  as 
the  independent  variables  are  replaced  by  the  velocity  magnitude  and 
direction  as  the  independent  variables,  linearizes  the  governing 
equation  which  may  then  be  analytically  solved.  The  plane  where  the 
said  velocity  parameters  are  the  independent  variables  is  termed  the 
hodograph  plane.  Likewise,  the  transformed  differential  equation  is 
known  as  the  hodograph  equation  and  it  is  solved  by  hodograph  methods. 

This  mathematical  study  addresses  the  solution  of  transonic  flow 
phenomena  by  an  extension  of  Lighthill's  hodograph  method. 
Lighthill's  method  transforms  a  given  solution  of  the  Laplace  equation 
into  a  solution  of  the  hodograph  equation  for  subsonic  flows  only.  A 
new  relation  is  developed  in  this  study  extending  this  transformation 
technique  to  include  flows  up  to  Mach  2.2735  in  air.  Requiring  only 
numerical  data  concerning  the  velocity  field,  this  hodograph  method  is 
computationally  efficient  and  mathematically  straightforward. 
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Chapter  1 


INTRODUCTION 


Problem  Overview 

The  differential  equation  in  terms  of  the  potential 
function,  <p  ,  governing  the  steady,  two-dimensional,  inviscid, 
irrotational,  isentropic  flow  of  an  ideal  gas  is 
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where  the  local  speed  of  sound,  c  ,  is 
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In  terms  of  the  stream  function,  4>  ,  (1.1)  and  (1.2)  are 
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d.n: 


For  a  derivation  of  the  above  results,  consult  Appendix  A 
and/or  Shapiro1.  The  non-linear  behavior  of  differential  equations 
(1.1)  and  (1.3)  is  evident. 

In  the  limit  of  vanishing  compressibility,  c2 — ) 


and 


equations  (1.1)  and  (1.3)  reduce  to  the  Lap-lace  equations 

^xx  +  ^yy  =  0 
^xx  +  ^yy  =  0  • 

The  potential  transform  of  Legendre  and  the  Molenbroek- 
Chaplygin  transformataion  linearize  equations  (1.1)  and  (1.3).  This 
study  addresses  the  solution  of  the  partial  differential  equation 


resulting  from  the  latter  transformation. 

Considering  the  physical  Cartesian  coordinates  x  and  y 
as  functions  of  the  potential  and  stream  functions,  one  may  deduce  the 


hodograph  equations 


=  q(pQ/p)  ;  =  q(p0/pq)<J>6 


(1.5) 


Defining, 


t  a  (q/qm)2 


where  qm  is  the  maximum  velocity  attainable  when  all  the  flow's 
internal  energy  is  converted  to  kinetic  energy,  the  governing 


hodograph  equation  is 


PQH>TT+PQT*T“*ee  =  0 

PQ  =  4 T2( 1 -t)/ (  1  ~t/ Tg) 


PQ  =  -Ht[1+(2-Y)/(Y-1)]/(1-t/ts) 
ts  s  (Y-1)/(Y+1)  =  Mach  one. 

The  above  equations  are  derived  in  Appendix  B. 


(1.6) 


(1.8) 


(1.9) 


The  solution  of  equation  (1.6)  presents  two  distinct 
difficulties.  First,  (1.6)  is  a  partial  differential  of  the  mixed 
type.  Consider  its  discriminant: 

B2  -  UAC  -  0  -1»(PQ)(-1)  =  HPQ  . 
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In  subsonic  flow,  0<t<ts  ,  the  discriminant  is  negative  and  hence 

equation  (1.6)  is  elliptic.  In  supersonic  flow,  ts<t<1  ,  the 

discriminant  is  positive  and  (1.6)  is  now  hyperbolic.  This  change  in 

behavior  elaborates  the  solution  method  considerably.  Second,  the 

Jacobian  of  the  Mol enbroek-Chapl ygi n  transf ormation  may  be  zero  or 

infinite  at  specific  points  or  along  lines  in  the  hodograph  plane. 

The  transformation  is,  therefore,  no  longer  one-to-one.  See  Appendix 

B  and/or  Ferrari  and  Tricomi2  and  Lighthill^  for  further  remarks 

concerning  the  Jacobian  of  the  transformation. 

Chaplygin,  in  1904,  solved  equation  (1.6).  The  solution,  a 

function  of  t,  0  ,  and  arbitrary  complex  separation  constant,  n  ,  is 

if>n(T,e)=Tn/2F(an,bn;n+1  ;i)e±ine  .  (1.10) 

The  hypergeometric  function  is  denoted  as  F(an,bn; n+ 1 ; t )  and  the 

product  xn/2F(an,bn; n+1 ;t)  is  termed  the  Chaplygin  function.  The 

reader  is  referred  to  Appendix  B  for  the  derivations  and  details  of 

(1.10).  Bergman*1  noted  that  any  particular  solution  (1.10) 

corresponding  to  a  particular  n  will  generally  converge  in  only 

part  of  the  hodograph  plane  and  that  the  complete  solution  for  the 

entire  hodograph  plane  is 

00 

<K  t  ,  e)=im{  £ipn(  Tce) )  (i.ii) 

n=o 

Appendix  C  and/or  Boerstoel-’  should  be  consulted  for  details  of  the 
hodograph  plane. 


In  light  of  Bergman's  result  and  of  the  fact  that  all 
boundary  conditions  are  lost  under  the  Mol  enbroek-Chapl  ygi  n 


transformation,  Lighthill^  posed 

CO 

4>(  T0)=lm{  Zi(/n(T)  f  (n,  Tm)  e-in®  (1.12) 

n=o 

as  a  solution  to  equation  (1.6).  The  normalizing  function,  fCn.t^), 
is  chosen  such  that  (1.12)  tends  towards  a  solution  of  the  Laplace 
equation  in  the  limit  of  vanishing  compressibility.  This  Laplace 
solution,  in  effect,  represents  the  boundary  conditions.  Lighthill 
developed  a  method  where  a  solution  to  the  Laplace  equation,  given  in 
terms  of  the  hodograph  variables  q  and  9  ,  is  transformed  into  a 
solution  of  the  hodograph  equation  (1.6).  For  subsonic  flow,  the 
transformation  is  especially  elegant,  requiring  only  numerical  data 
concerning  the  incompressible  velocity  field.  Supersonic  flow  regions 
require  a  Laurent  series  expansion  representation  of  the  Laplace 
solutions  in  terms  of  q  and  0  .  Nieuwland^  generalized  the 
representation  to  include  Mel  1  in-Barnes  integrals.  Both 
representations  require  highly  advanced  mathematics  to  properly 
represent  only  the  simplest  of  Laplace  solutions  in  the  hodograph 
plane.  In  fact,  the  Laplace  solution  governing  the  incompressible 
flow  about  a  general  lifting  airfoil  cannot  be  represented  by  either 
method.  The  trai 1 ing-edge  closure  problems  experienced  by  all  who 
employ  the  hodograph  method  are  symptomatic  of  this  fact. 

The  key  to  this  dichotomous  behavior  of  the  transformation 
lies  in  Lighthill’s  asymptotic  forms  of  the  Chaplygin  functions  as 
[ n j — *®  .  Because  equation  (1.6)  is  elliptic  for  0<t<ts  and 
hyperbolic  for  t3<t<1  ,  different  asymptotic  formulae  are  required  in 


each  region.  While  the  assumptions  behind  Lighthill's  subsonic 
asymptotic  formula  are  physically  sensible,  those  behind  his 
supersonic  asymptotic  formulae  are  not. 

Objective 

A  physically  reasonable  asymptotic  form  of  the  Chaplygin 
functions  as  |n| — and  valid  for  (Y-1  )/(Y+1  )<t<.5083  in  air  is 
developed.  A  transformation,  entirely  analogous  to  Lighthi 1 l's  is 
derived  where  only  numerical  data  concerning  the  incompressible 
velocity  field  are  required.  This  permits  the  transformation  of  far 
more  complicated  and  physically  realistic  Laplace  solutions  (hereafter 
termed  "model  flows"),  than  was  previously  possible.  In  particular, 
incompressible  flows  about  closed  lifting  profiles  can  be  transformed 
to  represent  transonic  compressible  flows  about  affinely  related 
profiles.  Lastly,  a  numerical  method  is  developed  which  calculates 
the  required  data  about  any  given  profile. 


Chapter  2 


AV. 


LITERATURE  SURVEY 

In  1890,  the  Dutch  mathematician  Molenbroek  linearized  the 
governing  equations  of  motion  by  considering  the  potential  and  stream 
functions  as  functions  of  the  velocity  coordinates. 

Chaplygin  in  1904  derived  the  solution  to  the  most  studied 
version  of  the  linearized  partial  differential  equation.  In  this 
method  the  strema  function  is  represented  by  an  infinite  series  of 
particular  solutions  each  of  which  converges  in  part  of  the  domain  of 
the  flow.  The  series  representing  the  strema  function  is  the  product 
of  a  hypergeometric  series  and  a  velocity  magnitude  parameter. 
Manipulation  of  such  functions  requires  rather  advanced  mathematics. 
Chaplygin  noted  that  by  specifying 


(pQ/p)  \/ M2-1  =  1 

the  equation  of  motion  reduces  to  the  Laplace  equation.  This  amounts 
to  specifying  a  fictitious  gas  with  Y  =  -  1  . 

Meyer  in  1908  found  a  "lost  solution"  which  is  the  expansion 
of  a  flow  around  a  corner. 

Demtchenko  and  Busemann  in  1932  and  1937,  respectively, 
realized  that  Chaplygin's  suggestion  of  Y  =  -  1  amounts  to  replacing 


7 


the  curve  of  the  ideal  gas  isentrope  in  the  p  vs.  1  / p  plane  by  a 
straight  line.  They  limited  themselves  to  considering  this  straight 
line  as  a  tangent  to  the  true  isentrope  at  stagnation  conditions. 
This  restricted  their  work  to  flows  of  Mach  number  less  than  about  .3. 

Von  Karman  and  Tsien  extended  the  above  method  by  taking  the 
tangency  point  to  coincide  with  free  stream  conditions.  This  shifted 
the  range  of  usefulness  to  higher  subsonic  Mach  numbers. 

Tsien  in  1939  appears  to  be  the  first  to  consider 
transforming  solutions  of  Laplace's  equation  to  solutions  of  the  mixed 
equation.  He  derived  a  relation  between  a  line  element  in  the  complex 
physical  plane  and  the  compressible  complex  potential.  His  last  step 
is  the  replacement  of  the  compressible  complex  potential  by  that 
corresponding  to  the  incompressible  case.  Thus,  a  transformation  from 
the  incompressible  physical  plane  to  the  compressible  plane  via  the 
hodograph  was  achieved. 

Ringleb  in  19*11  found  another  which  is  the  flow  about  a 
sharp  edge. 

In  1 9*4 6  and  19*17,  Tomotika  and  Tamada^  developed  an 
approximation  to  the  ideal  gas  isentrope.  Their  "isentrope"  coincides 
with  the  ideal  gas  isentrope  precisely  at  the  sonic  point  and  to  the 
order  of  their  tangents  at  the  stagnation  point.  The  resulting 
partial  differential  equation  for  the  stream  function  in  the  hodograph 
is  easily  solved.  The  solution  is  the  product  of  Bessel  functions  and 
trigonometric  functions.  They  calculated  a  shock-free  transonic  flow 
about  a  certain  airfoil-like  obstacle  without  circulation. 


8 


V"  \  ■- 


Nocilla  in  195^-6  used  the  Tomotika-Tamada  gas  to  calculate 
the  shock-free  transonic  flow  about  airfoils  with  blunt  noses. 
Circulation  is  again  absent. 

There  are  numerous  examples  of  researchers  approximating  the 
true  ideal  gas  isentrope  by  another  curve  rendering  the  solution 
easier  to  manipulate  than  the  Chaplygin  functions.  One  of  them,  the 
Tricomi  equation  replaces  the  ideal  gas  isentrope  with  a  straight  line 
and  the  regular  solutions  are  Airy  functions. 

Bergman^  in  19^5  appears  to  be  the  first  to  successfully 
treat  the  exact  case.  The  method  is  based  on  linear  integral 
operations  and  is  quite  complex.  Furthermore,  circulation  is  still 
not  included. 

Cherry®  in  19^7  solved  the  exact  equation  resulting  not  from 
the  Molenbroek-Chaplygin  transformation  but  from  the  Legendre 
potential  transformation.  In  this  transformat  ion  correspondence 
between  the  physical  plane  and  the  hodograph  ismore  direct  than  in  the 
Molenbroek-Chaplygin  transformation,  but  the  unknown  quantities  lack 
the  physical  basis  of  the  stream  and  potential  functions  present  in 
the  Molenbroek-Chaplygin  transformation.  Cherry  calculated  the  non- 
circulatory  flow  about  a  circular  cylinder. 

Lighthill®  in  19^7  is  the  first  to  solve  the  exact  equation 
resulting  from  the  Molenbroek-Chaplygin  transformation.  He  first 
developed  the  asymptotic  formulae  for  the  Chaplygin  functions  as  the 
magnitude  of  the  complex  separation  constant  tends  towards  infinity. 
For  subsonic  flow  only,  Lighthill  built  an  entire  function  in  the 


complex  plane  of  the  separation  constant.  This  entire  function  is 
composed  of  two  terms  which  turn  out  to  be  equal  to  one  another  by 
consideration  of  poles,  residues,  the  maximum  modulus  theorem. 
Lastly,  Lighthill  developed  the  compressible  stream  function  through 
a  transformat  ion  of  the  incompressible  complex  potential. 

That  this  method  even  exists  is  due  to  the  form  and 
simplicity  of  the  subsonic  asymptotic  formula  for  the  Chaplygin 
functions.  Lighthill  showed  that  if  a  general  Laurent  series  for  the 
Incompressible  flow  exists  in  the  hodograph  and  is  convergent  at  the 
sonic  speed,  then  the  equation  for  the  compressible  stream  function  in 
subsonic  flow  may  be  analytically  continued  into  the  supersonic  region 
of  the  hodograph  plane.  Such  a  Laurent  series  representing  the 
Incompressible  complex  potential  about  a  practical  airfoil  in  the 
hodograph  plane  will  be  very  difficult  if  not  impossible  to  derive  and 
manipulate. 

In  L ighthi 1 1' s  transformat  ion  techni que,  a  normalizing 
function  is  introduced  that  forces  the  solution  of  the  mixed  governing 
equation  to  reduce  to  solutions  of  Laplace's  equation  in  the  limit  of 
vanishing  compressibility.  Lighthill's  method  can  treat  flows  with 
circulation  via  the  requirement  that  as  one  encircles  the  airfoil  in 
the  hodograph  plane,  all  transformed  variables  remain  single- valued. 
This  limits  the  choice  of  the  normalizing  function  to  one  particular 
example  which  Lighthil 1  found. 

Nieuwland^  in  1  967  showed  that  Mel  1  in-Barnes  integral 
representations  of  the  incompressible  complex  potential  are  valid,  and 


preferred,  alternatives  to  the  Laurent  series  representations. 
Nieuwland  employed  the  Mel  1  in-Barnes  integrals  to  properly  model 
considerably  more  complicated  flows  than  could  be  handed  by  a  Laurent 
series.  He  calculated  the  flow  about  a  family  of  quas i- el  1 i pt i cal 
airfoils  with  circulation. 

Bauer,  Gorabedian,  and  Korn^  in  1972  solved  the  transonic 
flow  problem  by  rewriting  the  two  basic  governing  equations  in  a 
complex  form.  These  two  equations  are  then  decoupled  into  a  system  of 
linear  ordinary  differential  equations  which  they  solve  by  finite 
differences.  Bauer  et  al.  show  they  can  compute  a  wide  variety  of 
advanced  airfoils. 

Boerstoel^  in  1977  further  extended  the  method  of  Nieuwland. 
He  used  the  incompressible  complex  potential  Nieuwland  employed  to 
generate  a  "basic  stream  function"  which  possessed  the  required 
singularities  and  basic  properties  of  the  incompressible  complex 
potential  about  any  general  airfoil.  To  the  basic  stream  function, 
Boerstoel  added  the  "additional  stream  function"  which  also  satisfied 
the  governing  equation  and  represented  the  complex  incompressible 
potential  about  a  non-lifting  circular  cylinder.  An  assumed  arbitrary 
airfoil  image  in  the  hodograph  plane  permits  solution  of  the 
additional  stream  function  by  a  Tricomi  boundary  value  problem 
treatment  such  that  the  sum  of  the  basic  and  additional  stream 
functions  yields  the  pre-chosen  airfoil  image.  Boerstoel's  method 
produces  realistic  airfoils. 


correctly  reasoned  that  the  problem  lies  in  the  model  flow,  which  is 
transformed.  The  "YC"  profile  does  not  solve  the  problem.  It  is, 
however,  an  extremely  valuable  advance.  With  relative  ease,  the 


transonic  flow  about  a  wide  variety  of  cambered  lifting  airfoils  is 


calculated.  These  airfoils  are  similar  to  Joukowski  airfoils. 


Chapter  3 

ASYMPTOTIC  FORMULA  DEVELOPMENT 


Lighthill's  incompressible  flow  to  subsonic  compressible 
flow  transformation  depends  upon  two  facts.  First,  any  solution  of 
the  hodograph  equation  (1.6)  must,  in  the  limit  of  vanishing 
compressibility,  reduce  to  a  solution  of  the  Laplace  equation  in 
hodograph  variables.  Second,  the  maximum  value  of  a  function  analytic 
and  not  constant  in  a  domain  occurs  on  the  boundary  of  that  domain  by 
the  maximum  modulus  theorem. 

Consider  the  first  point.  For  strictly  subsonic  flow,  the 
characteristic  equation  (see  Appendix  D)  is 

d9/aT=(-PQ)_1/2  (3.1) 


which  when  integrated  yields: 

s  =  o  tanh"1  -iL+tanh-1  /  W^) I 

Vy-1  Y  (y+1)  (1-t)  \J  (y-D(l-r)  (3.2) 

where  0  is  an  arbitrary  constant.  The  hodograph  equation  for 
strictly  subsonic  flow  is 


^ss+,J,09=T(s)^ 


(3.3) 


wh©r6 

T(s)=  -y-(-PQ) 1/2 (TT  -  — )=  t2(1-t)"1/2(1-  ^)3/2 

2  Q  P  (y-1)  y_1  (3.M 
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Thus,  s  and  T  are  both  analytic  functions  of  t  in  the  domain 

t < Tg  .  Since  the  derivative  of  s  with  respect  to  t  is  not  zero  in 

the  domain,  t  is  an  analytic  function  of  s  by  the  inversion 

theorem.  Since  x*xs  or  1  in  said  domain,  T  is  also  analytic 

function  of  s  .  For  very  small  x  ,  the  Chaplygin  functions  behave 
as  •  Because  the  Chaplygin  functions  comprise  at  least  part  of 

the  compressible  flow  solution,  the  functional  form  of  x  with 
respect  to  s  is  needed  so  that  the  compressible  solution  can  be 
forced  to  the  Laplace  solution  in  the  limit  of  vanishing 
compressibility.  Lighthill^  chose  the  value  of  o  to  be  that  value 
which  causes  "t  to  be  asymptotically  e  as  t— >0  and  s — *•-»  ,  o 
is  that  value  of  s  at  x  =  (Y-1  )/(Y+1 ) ,  the  sonic  speed;  .  . 
That  value  of  a  is 

a=  j£n|  2(y-l)  tanh  1 

Details  of  the  derivation  are  in  Appendix  D. 

Lemma  1  (Lighthil  1^).  Both  t  and  T  are  analytic  functions  of  e23 
in  the  region  |e2sj<e2°  . 

p  o 

Proof  .  .  .  e  =t  ...  is  an  analytic  function  of  x  in  any  domain 
excluding  x=(Y-1 )/(Y+1 )  and  t=1  ;  its  derivative  is  not  zero  in  such 
a  domain;  and  at  neither  of  the  singular  points  can  je2s|  be  less 
than  e2°  .  Hence,  by  the  inversion  theorem,  x  is  an  analytic 
function  of  e2s  in  the  region  ]e23j<e2°  .  As  x  is  never  1  or  (Y- 
1 )/( Y+1 )  is  the  region,  T  must  also  be  analytic. 

For  strictly  supersonic  flow,  the  above  incompressible 


o=  -1.173  for  o=1.4. 


(3.5) 


boundary  condition  does  not  exist.  Integrating  the  characteristic 
equations  of  supersonic  flow  yields: 


t  =  +{eH 


Y+l  .  ,-1,  l(.+i)  ■-(•-!)  -1,  (i~H) 

tanh  \l  ,  . , N  ,,  x  \  (,-!)(!-:) 


Y-i  —  y  (y+d  (i-o 

while  the  hodograph  equation  for  purely  supersonic  flow  is 
*tt+S(t)^t=^60 


(3.6) 


(3.7) 


where 


s(t)  = 


-  (PQ) 


1/2  Or 
C  Q 


2  (  '  +  1 ) 
(.-l)2 


2  -1/2 
(1-  )  ' 


.  +  1 
-1 


:-i) 


-Ml 


(3.3) 


The  value  of  e  is  that  value  such  that  t  is  asymptotically  e±2l"1 
as  x — *(Y-1)/(Y+1)  .  At  a  glance 


e=1 /2£n | Y-1 /Y+1 |  ;  e=  -.8959  for  Y  =  1.4  .  (3.9) 

Appendix  D  contains  further  details  concerning  the  supersonic  flow 
equations  where  t  =  t^e  . 

+  9  4- 

Lemma  2.  Both  t  and  S  are  analytic  functions  e~  in  the  domain 
ts<t<1  . 

Proof.  By  equation  (3.6),  e±2t1  =  t+  .  .  .  is  an  analytic  function 

of  t  in  the  domain  ts<t<1  .  No  singular  points  exist  in  this 

domain  and  its  derivative  is  not  zero  in  the  domain.  By  the 

inversion  theorem  x  is  an  analytic  function  of  e±2t1  in  said 

domain.  Because  e2e  is  a  constant,  then  t  is  also  an  analytic 
+  ?  f 

function  of  e-  .  Since  S  is  an  analytic  function  of  x  ,  it  is 
also  an  analytic  function  of  e±2t  . 


Consider  now  the  implications  and  requirements  of  the 


maximum  modulus  theorem.  The  Chaplygin  functions  are  analytic  with 
respect  to  n  except  at  n  =  -2,  -3,  where  they  have  simple 

poles . 

Theorem  1  (Lighthi  1 1^).  If  0<xS1,  vpn( t )  is  an  analytic  function  of 
n  except  at  n  =  -2,  -3,  -4,  .  .  .  ,  where  it  has  simple  poles,  its 
residue  at  n  =  -m  being  ,  where  Cm  .  .  .  is  positive  and 

~  (2irm)_1  e""2om  as  m — ►»  . 

The  sign,  ~  ,  means  "asymptotes  to."  Assume  two  functions  defined  on 
the  entire  complex  plane  can  be  constructed  which  possess  identical 
poles  and  residues  at  those  poles.  Their  difference  will  be  an  entire 
function.  The  maximum  value  occurs  at  infinity  by  the  maximum  modulus 
theorem.  If  this  maximum  value  is  zero,  by  the  maximum  and  minimum 
modulus  theorems,  the  two  functions  are  equal  to  each  other  on  the 
entire  complex  plane.  Lighthill  employed  this  fact  in  his 
incompressible-to-compressi  bl  e  transf  ormation  where  one  of  the  two 
functions  is  constructed  of  Chaplygin  functions.  The  forms  of  the 
Chaplygin  functions  as  |n| — ,  called  asymptotic  forms,  are 
required  to  ensure  the  maximum  modulus  is  zero.  The  asymptotic  forms 
change  as  the  flow  changes  from  subsonic  to  supersonic. 


Subsonic 


Integrate  T(s)  to  produce  V(x) 


Thus, 


V(t)=(~P/Q) '  =EXP{1 /2 


V(0)=1  ;  dV( t )/ds=TV/2 


)ds,}  . 


(3.10) 


(3.H) 


Assuming 


!^n(T)=ensV(  f)Wn(s)  (3-12) 

the  hodograph  equation  for  strictly  subsonic  flow  (3.3)  becomes 

d2Wn/dS2+2ndWn/ds=[1/4T2(s)-1/2dT(s)/ds]Wn  .  (3.13) 

Lighthill  determined  that  as  | n | ■ — and  excluding  the  negative 
integers,  W  (s) — *1  ■ 

Theorem  2  ( Lighthi  1 1  ^).  If  6>0  and  o1  <o  ,  then  Wn(s)  1  ,  i.e. 
i|)n(T)~ensV(T)  ,  uniformly  for  s^o^  and  for  n  in  the  whole  complex 
plane  with  circles  of  radius  6  around  each  negative  integer 
excluded,  as  |n| — ►<»  . 

Consult  Appendix  D  for  the  complete  derivation  of  theorem  2. 

Supersonic 


where 


For  ts<t<1  ,  define: 

fl(  t)  =  (P/Q)1A+  constant=EXP {  -1  /2  /s(t2)dt2) 


(3.14) 


fl(xs)-1  ;  dn(T)/dt=  -Sn/2  .  (3.15! 

Assume, 

<i/n(T)=entfl(T)  Ln(t)ene  .  (3.16) 

Since  <|>n(T)  ,  ent  ,  and  fl(x)  are  analytic  functions  of  t  and  e2t 

by  lemma  2,  so  too  must  Ln(t)  .  Thus, 

00 

Ln(t)  =  I  ln,re2rt  *  (3.17) 

r=o 

Note  from  equation  (3.16)  lim  Ln(t)ene  =  «|>n(ts)  •  Subjecting  equation 

t-»o 

(3.16)  to  the  linear  operator  of  equation  (3.7)  yields: 
d2Ln(t)/dt2  +  2ndLn(t)/dt  = 


(3.18) 


[1  /4S2(t)+1 /2  dS(t)/dt  -2n2]  Ln(t) 

By  lemma  2 

00 

1/4S2(t)+1/2  dS( t )/dt  =  Z  Spe2rt  .  (3.19) 

r=o 

Substituting  equations  (3.17)  and  (3.19)  into  (3.18)  yields: 


,  2,  2rt  „  „  .  2rt  _  2rt  2 

L  4r  e  +2n  Z  2r l  e  =  Z  s  e  -2n  < 

,  n,r  .  n,r  r  t 

r=l  r=l  r=o  r=o  n,r 


2rt 


(3.20) 


Pt 

Equating  powers  of  e 

4r (n+r) i 


L  S  l 


■zb. 


n’r  m=o  r-m  n,m  ~‘n,r 


(3.21  ) 

2 

-  "  s  i  +  (s  -2n) i  .  (3.22) 

m=o  r-m  n,m  o  ’  n.r 

Equating  the  e°  term  yields  SQ  =  2n2  or  ln  Q  =  0  .  Assuming  £n  Q 
=  0  implies  that  any  general  ln  p  is  a  function  of  1R  1  .  Since 
ln>1  is  not  restricted  in  any  way,  S(t)  is  implied  to  be  a  function 


of  an  arbitrary  constant  which  is  false.  Thus,  SQ  =  2n"-  and 


Hr (n+r )l 


r-1 

=  r  s  i 
n,r  Jr-mvn,ra 
m=o 


(3.23) 


This  equation  predicts  the  occurrence  of  poles  at  negative  integral  n 
for  (iJn(f)  as  Hn  r  becomes  indeterminate  at  these  points.  At  first 
glance,  the  assumption  of  SQ  =  2n  would  seem  to  imply  S(t)  =  f(n) 
which  is  false.  However,  the  implication  that  the  coefficients  Sp 
are  functions  of  n  does  not  guarantee  that  S(t)  is  a  function  of 
n.  Lighthill  also  must  use  this  observation  in  his  development  of  the 
asymptotic  formula  of  the  Chaplygin  functions  for  subsonic  flow.  By 
lim  Ln(t)eni=  ij>n(Ts)  equation  (3-16)  predicts 


(3.25) 


Ln(t)eIIt'  <  ^n(ts)  E  e 
r=o 

By  the  triangle  inequality, 

Ln(t)ene  -  iJ»n(Ts)  £  iJ>n(Tg)  E  e2rt  .  (3.26) 

r=1 


Assuming 

♦n(T)  -  e~nt  n(T)Hn(t)ene  (3.27) 

results  in  equations  similar  to  equations  (3*27)  through  (3.26)  with 
t  substituted  by  -t  throughout.  Thus,  in  general, 

*n(T)  -  «(T)ene[Ln(t)ent+Hn(t)e"nt]  (3.28) 

! Ln( t )enE |  -  |«n(Tg)|  S  |*n(Tg)  E  e2rt| 

r=1 

00 

|Hn(t)ene|  -  |*n(Tg)|  S  |«Pn(Tg)  E  e'2rt|  .  (3.29) 

r=1 

Equations  (3.26)  and  (3.29)  can  be  investigated  more  thoroughly  with 
the  following  theorem. 

Theorem  3  ( Lighthill^) .  When  |arg  n|  £  m  -  5  , 

<J»n(Y-1  /Y+1 )  ~  kenon1 

where  k  =  it  B(  1  /3[T( 2/3 )  1_1 ) 2( 2Y-1  )/(2Y-2)  ( Y+1  )-(  Y+2)/(6y-6)  . 

When  |arg  ( —  n )  |  ^-rr— 6  ,  however, 

♦n(Y~1/Y+1)  ~  keno(-n)1/^  sin(mr-ir/6)/sin(mr)  . 

Equations  (3.26)  and  (3.29)  hold  only  when  I 2t I  £  I  o  I  ,  otherwise  the 


right  hand  sides  are  greater  than  £n(-rs).  Thus  for  |t|<|o/2|  , 

eneLn(t)  =  eneHn(t)-^^n(Ts)  .  (3-30) 

Theorem  4.  For  |t|<|o/2|  and  excluding  circles  of  radius  6  about 
the  negative  integers, 

ipn(x)  ~  keno|n|1/6fi(T)[ent+e~nt]X(n) 
where  Xn=1  when  |arg  nj£  it —6 

Xn=  sin(mr-ir/6)/sin(niT)  when  |arg  ( —  n )  |  ^  tt — 6  as  |  n  |-“  . 
The  stipulation  that  |t|£|o/2|  =  .5867  corresponds  to 


t  =  .5083  and  to  a  maximum  Mach  number  of  2.2735  for  Y  =  1.4. 


Chapter  4 


THE  INCOMPRESSIBLE-TO-COMPRESSIBLE  TRANSFORMATION 


Transformation  in  the  Hodograph  P lane 

The  series  representing  each  x )  converges  only  in  part 

il  o 

of  the  hodograph  plane.  Bergman  ,  Bers  and  Gelbart,  and  Lighthill-3 
thus  assume  the  total  compressible  stream  function  can  be  constructed 


from  an  infinite  series  of  ^n(t)  , 

CO 

I  <l>n(T)e”in  9}' 


(4.1) 


Consider  an  incompressible  complex  potential 


1  ^  ^  and  define 

JA(i) 


u-iv  -  qe 


"iO  -  s 


(4.2) 


(see  Appendix  C).  If  such  a  complex  potential  is  analytic  outside  the 


body,  then  in  hodograph  variables  *>  *> 

(i)  \  '  n  \  '  n  -ii 

’  ‘  2^  V  ’  Zj  Cnq  e 


♦<l)  -  Ira  (  ^ 


n  -  inO 
c  q  e 
n 


(4.3) 


(4.4) 


The  similarity  of  equations  (4.1)  and  (4.3)  prompted  Lighthill  to 


assume 


E 


*  *  Im  {  2^Vn(T)f<n’\,>e  >  (4.5) 

n=o 

where  fCn.T^)  is  called  the  normalizing  function.  Its  primary  duty 
is  to  force  equation  (4.4)  to  show  the  same  behavior  as  equation  (4.5) 
in  the  limit  of  vanishing  compressibility. 

In  particular,  for  o  <  Too  <  ts  Lighthill  stipulated: 

(a)  fCn.Ta,)  is  an  analytic  function  of  n  except  possibly  at 
certain  real  non-negative  poles  of  each  of  which  it  has  a 
real  residue. 

(b)  For  large  n  fCn.t,,,)  ~  Ae  ns“  for  some  constant  A  . 

(c)  As  — *o  ,  fCn.x*,)  ~  Tm_n/2  uniformly  for  all  n  at  a 

distance  6  from  any  pole  of  fCn.r^)  . 


Condition  (b)  ensures  parallel  behavior  between  equations  (4.5)  and 
(4.4)  for  large  n  and  t  =  t.  .  Equation  (4.5)  for  large  n  and 
o  <  t  <  ts  becomes:  00 


i|»  =  Im  { 

T  c  V(x  )e_inf'  } 
x  j  n  °°  1 

(4.6) 

n=o 

while  equation  (4.4)  with 

q  =  1  is  <» 

..(i)  T  r  \  >  -inO 

V  =Im  t  /  cne  } 

(4.7) 

n=o 

Condition  (c)  ensures  that 

equation  (4,5)  reduces  to  equation 

(4.4)  as 

dpi-*”  .  In  this  case 

<J»n— i *Tn/2[1+0(x)] 

and  t|/n(-r)f  (n,  t.)  becomes  (t/On/2  =  (q/qj0  =  q11  for  qm 
normalized  to  one. 

Lighthill  considered  closed  contours  about  the  origin  of 
the  complex  n  plane.  For  t  <  ts  ,  Lighthill  postulated  a  function: 
J1  (n.T.x J  =  1  /nUn(T)f(n,TM)en(3»“s)-f(o,TM))  . 


(4.9) 


One  notes: 


(i)  Because  ipQ(  t  )  f  (o,TC0)e°=  f  (o,  rm)  a  removable 
singularity  exists  at  n  =  o  • 

(ii)  By  theorem  2,  for  large  n  ,  J1  ~  1/n(AV(t)) — »o  as 

I n |  • 

(iii)  By  condition  (a)  and  theorem  1,  the  poles  P1.P2.P3>  • 

of  i|/n(T)f(n,Tai)  are  real  as  are  all  the  residues 

r1  'r2,r3  •  •  •  of  n_1  *n(T)f(n,Ta>)  . 

Lighthill  postulated  another  function: 

00  ^  „  m(  300-s) 

J2(n>T»Too)  = 

One  notes: 

00 

(iv)  If  £  rmepn/s<»~3)  converges  absolutely  then  J2  is 

m=1 

an  analytic  function  of  n  except  at  the  poles 
Pi «P2»P3»  •  •  • 

(v)  As  | n | — *<*  and  n  remains  a  distance  6  away  from 
each  pm  ,  J2 — *0  . 

(vi)  At  each  pole  pm  the  principal  part  of  J2  is  the  same 
as  that  of  J1  . 

Lighthill  noted  that  the  function  J1  -J2  is  entire  in  the  n  plane 
by  observations  (i),  (iii),  (iv),  (vi).  By  observations  (ii)  and  (v), 
its  maximum  modulus  approaches  0  as  [n| — . 

Thus:  .  P„<s„-s> 

11  -  via! _ 


fu.m 


[(T)f(n,iB)  =  f  (o)Tjj)en(s'3°"s  }  +  n  ^  - 


(n-p_)  (s  -s  ) 


ra 


m=l 


n"P„ 


(4.12) 


and,  as  n| — *a>  yet  remaining  a  distance  5  from  each  pm  , 


AV(t)  =  f(o,Tj  +  ^ 


Pjs~  s).: 


r  e 
m 


m 


m~l 


(M. 1 3) 


Substitution  of  equation  (4.12)  into  equation  (4.5)  yields 


'Wm  f (o,tJ 


n  (3  -s  -i~)+ 


^  rrae'1Pm'‘  X) 


1 


m=l 


n-Pm  | 

m  J 


(4.14) 

The  form  of  the  left-most  summation  suggest  £  =  e3-s“-10  when 
compared  with  equation  (4.4).  Thus 


s  —  s„— 


it; 


oo 

E  / 


i|i=Im  •{  f(o,Tj^^  +  y  r  e  lpm"  /  nc  pm~1dr 
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noting  d<J>^  =  Enc  5n  ^dc  from  equation  (4.3)  'o 


Lighthi  1 1  noted  that  gm  =  f(m,T00,co)  •  the  pm  is  a  pole  of 
4<n  ,  then  rm  =  f(i|>_pm)  by  theorem  1.  If  the  pm  is  a  pole  of 


f(n,T,„)  ,  then  rm  =  r(<4pm)  •  Thus 
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is  independently  a  solution  of  the  equation  of  motion  as  is  the 


remainder  of  equation  (4.16),  or 
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Assuming  pQ  =  o  ,  rQ  =  f(o, ;  )  ,  and  taking  the  imaginary  part  of 
equation  (4.18),  yields: 
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Pm  (-s inp  (  —  )d-tl)+cosp  (  -•  )(1  (i) 
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The  subscript  one  denotes  the  variable  of  integration  and  £0  »  Q0  *  0 
for  pm  >  1  to  ensure  convergence  of  equation  (4.18).  The  relation 
z,  =  es-s<x>_19  is  the  transformation  from  the  incompressible  flow  plane 
to  the  compressible  flow  plane.  Thus  q  =  es_s°°  defines  the  speed 
magnitude  transformation  while  the  flow  angle  0  is  invariant  under 
the  transformation. 

Equation  (4.14)  transforms  incompressible  flow  solutions  to 
compressible  flow  solutions  via  simple  numerical  integration  if  the 
details  of  the  governing  Laplace  solutions  are  known  numerically. 

An  equation  similar  to  equation  (4.19)  is  now  derived  for 
supersonic  flow.  The  procedure  is  exactly  analogous  to  Lighthi  1 l's 
subsonic  case.  Note  also,  no  constraints  exist  on  f(n,0  with 
respect  to  t  .  Thus,  any  form  of  f(n,ia;)  valid  for  subsonic  flow 
is  also  valid  for  supersonic  flow. 


where 


Consider,  for  0<|tj<|o/2|  ,  t  >  0  , 


J3(n,x,T (vn(T)f(n,gen(Va  t)-E(o,uiu)) 
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(4.20 


Again,  note: 

(vii)  A  singularity  of  order  1/6  exists  at  n  =  0  . 

(viii)  By  theorem  4  and  condition  (b)  for  t>  0  ,  |J,[ — *0  as 

| n | — while  remaining  a  distance  6  from  any  pole  pm  . 
(ix)  Observation  (iii)  holds. 

Consider,  for  o<|t|<|o/2|  ,  t>  o  , 


J/.(n>T»0  = 


sPm(soo-a-t) 


m  n-p 
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Note: 


A  singularity  of  order  1/6  exists  at  n  =  0  . 
r  and  p  each  include  an  additional  value  of  0  at  n  = 


(xiii)  Observations  (iv),  (v),  (vi)  hold. 

Thus  is  an  entire  function  in  the  complex  n  plane  whose 

modulus  tends  to  zero  as  jn| — while  avoiding  the  poles  pm  . 
Thus  =  Jjj  and  for  0<|t|<|o/2j  ,  t>  0 

,  .  "  (n-pm) (a+t-s„) 

*n(T)f(n,Tj=f(o,Tjen(°+t-s”)+ny'r  1 _ — _____ 
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and  as  |n| — *■«  while  remaining  a  distance  6  from  the  p 


|  n|  (x)  =  f  (0,t  )+  S  \ 
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(4.23! 
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Equations  (4.22)  and  (4.23)  -are  exactly  anologous  to 
equations  (4.12)  and  (4.13).  Thus,  by  inspection,  the  form  of 
equation  (4.18)  applicable  to  supersonic  flow  where  t  >  0  is, 


sIra 


f(o’:..);(i)+  X/V~iPm'  f  -’Pm  j  .  ( i ) 
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(4.24) 


Letting  pQ  =  0  ,  rQ  =  f(o,Tm)  and  taking  the  imaginary  part  of  the 


above  gives 
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where  one  denotes  a  variable  of  integration.  If  t  <  0  ,  then 
substitute  +t  into  equations  (4.20),  (4.21),  and  (4.23)  and  -t  into 
equations  (4.20)  and  (4.24).  However,  the  result  is  the  same  in 
either  case  as  -(-t)  =  +  t  .  Equation  (4.25)  is  therefore  taken  to  be 
the  supersonic  analog  to  equation  (4.19).  Now,  q  =  e0  +  t  s°°  and  0 
is  still  invariant  under  the  transformation. 

Passage  to  the  Physical  Plane 

Derivation  of  the  coordinates  of  the  physical  plane  x  +  iy 
=  z  begins  with  the  Molenbroek-Chaplygin  transformation, 
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(4.27) 


(4.28) 


P  /p=(1-t)  ’ 
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(4.39) 


Thus, 


p  /pQ  =  1/2t 
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(-4.30) 


Recalling 


Equation  (4.28)  becomes 
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(4.32) 


The  values  of  ij;g  and  (J»T  are  obtained  from  equations  (4.19)  and 
(4.25)  for  subsonic  flow  and  supersonic  flow,  respectively.  However, 
the  only  differeence  between  these  two  equations  lies  in  the 
relationship  between  q  and  t  .  Equations  (4.19)  and  (4.25)  yield 
exactly  the  same  relationships  when  operated  on  by  equation  (4.32). 

The  remainder  of  this  entire  chapter  is  from  Lighthill^  with 
the  exception  of  the  generations  to  include  supersonic  flows. 

oo 
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Redefine  V(r),  |n|  &fi(x)  as  £n(x)  to  Permit  the  development  of  a 

general  passage  to  the  physical  plane,  irrespective  of  flow  character. 


Note  that 


"  rm  q  =  A  £  (-)  by  equations  (4.13)  and  (4.23)  and  by 
the  incompressible-to-compressible  speed  magnitude  relations.  Next, 


use  the  facts 
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to  reduce  equation  (4.33)  to: 
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Noting  pQ/p  =  1  in  incompressible  flow, 
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Let  the  iucompressible-to-compressible  transformation  be  denoted  ek 
where  k  =  s-s^  subsonically  and  k  =  o+t-s^  supersonically. 
Equation  (4.36)  becomes, 
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The  first  line  of  equation  (4.40)  corresponds  to  the  first  and  third 
terms  of  equation  (4.36).  The  second  line  corresponds  to  the  second 
and  fourth  terms  while  the  third  line  corresponds  to  the  sixth  term 
and  the  last  line  corresponds  to  the  last  term. 

Collecting  the  coefficients  of  Z  yields  Q(t00/t)1/2  times 

1  a5  ek+  ek+  e  (r '+  J^) 

2  n  dx  4t  n  ^  m  2x  2(l-pm)  *  (4.41 

Redefine  equations  (4.12)  and  (4.22)  as 
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and  equations  (4.13)  and  (4.23)  as, 
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where  rQ  =  f(o ,xj  and  pQ  =  0  .  Equation  (4.42)  with  n  =  1  is: 
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Using  the  above  relations  with  equation  (4.43)  •  equation 


(4.41 )  reduces  to: 
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Similarly,  the  coefficient  of  £2dz  is  Q(t00/t) 


times 
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Equation  (4.47)  reduces  to 
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which,  by  the  way  of 


further  reduces  to 
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Finally, 
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The  arbitrary  function  of  t  ,  g(t),  is  found  to  reduce  to  an 
arbitrary  constant  upon  operation  under  the  counterpart  to  equation 
(4.32),  namely, 
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2r  r t 


v  ‘  (4.52 

This  exercise  is  not  followed.  The  reader  is  referred 


to  Boerstoel's3  concise  treatment  of  the  exercise. 

Equation  (4.51)  was  derived  under  the  assumptions  pm*  1  ,  - 
1.  By  theorem  1  and  condition  (a),  no  pole  exists  at  n  =  -  1  . 
However,  it  is  possible  for  f(n,TB)  to  have  a  pole  at  n  =  +  1  .  A 
new  equation  giving  the  correspondence  between  the  hodograph  plane  and 
physical  plane  and  valid  for  pm  =  1  must  be  deduced.  The  process  is 
very  similar  to  the  one  given  above  and  the  reader  is  referenced  to 
Lighthill^  or  Boerstoel'’  for  the  treatment. 
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Normalization  Functions 


The  actual  form(s)  of  f(n,T00)  must  be  found  satisfying  not 
only  conditions  (a),  (b),  and  (c)  but  also  ensuring  that 
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converges  absolutely.  An  obvious  choise  is 
f(n,T_)  =  e~ns°° 


(4.55) 


given  by  condition  (b).  Conditions  (a)  and  (c)  are  also  satisfied. 
Because  no  poles  exist,  the  convergence  of  equation  (4.54)  is 


By  theorem  1  , 


determined  by  the  rm  due  to  the  ^n(x)'s  . 
convergence  is  assured  if  (o  +  s^-sXO  for  subsonic  flow  and  if  (sM- 
t)<0  for  the  supersonic  case.  The  subsonic  case  is  always  true.  The 
supersonic  case  is  true  beyond  the  range  of  validity  of  theroem  4  if 
the  free  stream  is  subsonic.  Because  n  =  1  is  not  a  pole,  equation 
(4.51)  may  be  used  for  passage  back  to  the  physical  plane. 

Recall  that  the  hodograph  surface  possesses  a  singularity  at 
£  =  1  which  corresponds  to  the  free  stream  conditions.  Thus,  when 
circulation  is  present,  the  integrals  of  x,  to  any  power  with  the 
variable  of  integration  as  Z  will  all  suffer  a  fixed  increase  once 
this  free  stream  singularity  is  encircled. 

Derivation  of  the  normalization  function  for  circulatory 
flow  depends  on  forcing  the  said  integrals  to  remain  one  valued. 
Consider  equation  (4.16)  and  note 

d<t(i)  =  £dZ  .  (4.56) 

The  integral  in  equation  (4.16)  becomes  s1-Pm  dz  .  For  x,  near  one, 


C1  Pm=l+(l-p  )  (C-1)+0|  r,-l|  2 


(4.57) 


The  integral  ^^dz  is  zero  while  the  integral  /td  z  is  simply  the 
circulation,  r  .  Thus,  the  integral  ^1-Pm  dz  suffers  a  fixed 
increase  of  ( 1  —  pm) r  as  the  free  stream  singularity  is  encircled. 
Therefore,  equation  (4.16)  and  (4.25)  increased  by 
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This  sum  can  be  zero  only  if  the  pm's  are  symmetrically  placed  about 


the  origin  in  the  n  plane  (except  possibly  pm  =  1  which  needs  no 


counterpart).  Denoting  these  poles  as  p  ,  =  -p  ,  the  residues  must 
satisfy 
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The  residue  of  n  ^ij/n(x)f  (n,  x*)  due  to  i/in(x)  at  n  =  -m(for  m£2) 
is 

Cm^m^ T ) f  (_m>  "O  •  (4.60) 

By  equation  (4.59)  that  at  n  =  m  must  be 

( 1  +m/1  -m)CmiJ>m(  x)f  (~m,  t0)  (4.61) 

so  that  the  residue  of  f(n,T00)  at  n  =  m  is: 

(4.62) 

Lighthill  next  found  the  form  of  f(n,xa>)  which  satisfies  the 
above  as  well  as  conditions  (a),  (b),  and  (c)  through  educated 

guesses.  It  is: 
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A  pole  exists  at  n  =  1  .  Equation  (4.53)  should  also  be 
one  valued  with  this  form  of  f(n,x0B)  .  As  the  free  stream 


singularity  is  encircled,  equation  (4.53)  increase  by: 
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By  equation  (4.59)  the  summation  may  be  recast  as 
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which  is  zero  by  the  symmetry  about  n  =  0  .  Equation  (4.53)  is  one¬ 


valued  only  if 


By  equation  (4.63) 
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Substituting  all  of  the  above  into  equation  (4.66)  gives 
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Thus,  equation  (4.66)  does  indeed  equal  zero  and  one  can 


state  conclusively  that 


f (n,i  )= 
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is  a  proper  normalizing  function  for  circulatory  flows. 
At  the  poles: 


(i)  n  =  -m(m-2, 3, 4,...) ,  the  residue  of  n~  ip  ( x)f (n, j,)  is 


(4.74) 


(il) 

(iii) 


Lighthill 


n=m(m=2 , 3, 4 , . . . ) ,  the  residue  is 

n  =  1  ,  the  residue  is 


(  t)0_i  (  O  +2t0O^'  (tb)] 

3  stated  that  equation  (4.53)  becomes 

/" 


Z=F(t  ,  t  )Z+ 


T  (1-r  ) 


1/,-L 


>-I 


K-XdZ+i  Z- 


<iZ 


-1/t-I 


■>  i  •  ■  „  1  /  ? 

Zo  ]+Q(— ) 


•:  c  (.  (T  ) 

n  m  « 


(4.75, 


(4.76) 


^(n+l)  ~  1-n 

+2x^L(^ f-  .  (-))  b  ^  (  h - dZ 

*  ro  *  m  2-.  m  2(m+l)  J  , 

l-m 


/*n+l  i  ( 1-n)  •  /*  .-n+1 

dz)  +(,’(T)+  ~L,  (:))  - (  /  - — —  dZ 

ro+1  ra  2t  ra  2(l-m  J  m+1 


m+1 

L-cx 


-/w  dz)> 


(4.77) 


where 


F(r,Tj=i+-  a-a-xj 


y/y-1 


)+ 


1/y-l 


[i'ni 


Too  -Y/Y-l  T 

+a+i,  .(130  !_  -ldxt-L 


(4.73) 


o 

Finally,  by  theorem  1,  the  convergence  of  equation  (4.54)  is 
guaranteed  if  (2o-s+sM)<0  for  subsonic  flow.  This  is  never  violated. 
By  theorem  4,  convergence  is  secured  if  (o  +  s<JO-t)<0  which  is  again 
valid  beyond  the  range  of  validity  of  theorem  4.  Thus,  the  assumption 
of  convergence  required  in  the  transformation  is  valid  under  these  two 
choices  of  normalizing  functions. 

The  reader  is  again  referred  to  Boerstoel'*  who  derives 


Chapter  5 

CONSTRUCTION  OF  THE  MODEL  FLOW  AND  THE  VELOCITY  FIELD 

Construction  of  the  model  flow,  or  a  governing 
incompressible  complex  potential,  is  required  because  Lighthill's 
transformation  technique  requires  a  detailed  numerical  knowledge  of 
the  governing  imcompressible  complex  potential  and  the  associated 
velocity  field. 

The  coordinates  of  the  profile,  =  0  ,  to  be 

transformed  are  assumed  to  be  known.  Also,  it  is  assumed  that  a 
combination  of  doublets  and  clockwise  vortices  placed  in  a  uniform 
flow  represent  the  profile.  The  strengths  and  positions  of  the 
doublets  and  vortices  are  to  be  found  through  the  knowledge  of  the 
given  profile  geometry.  In  what  is  to  follow,  it  is  assumed  that  the 
doublet  centers  and  vortex  centers  are  by  no  means  coincident. 

Consider  the  complex  Z-plane  in  which  the  doublet  and 
vortex  centers  are  placed.  See  Figure  8.  Assume  a  Zj  coordinate 
center  is  placed  at  the  center  of  the  jth  singularity  -  be  it  a 
doublet  or  a  vortex.  Then,  at  a  general  point  P  ,  the  incompressible 
complex  potential  is: 
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+ 1  ^  -  q  z+  l  — —  +  —  I  1  ,  ,  n  ,  z 

,  Z  2'  ,  ,  b  b 

a=l  a  b=l 


(5.1  ) 


Setting  qm  =  1  and  normalizing  the  doublet  and  vortex  strengths 


yields: 


MK  (1>  =  7  +  7  -§■  +  -i 


Z  +  E  t  +  J^  E  7b;n!zb: 

a=l  a  b=l 


(5.2) 


Note  that  the  jth  singularity  center  coordinates  are  Zj  =  Rjei0j  in  the 
Z-plane.  The  coordinates  in  the  jth  Z-plane  are  Zj  =  rje*0^  .  Taking 
the  real  and  imaginary  parts  of  the  incompressible  complex  potential  gives: 


/  • \  n  B  cos i  ,  m 

U)  _  „  .  _  a  a  1 

I  =  RcosO  +  z  -  -  —  z  r.ot. 

,  r  2x  ,  _  b  b 

a=l  a  b=l 


=  RsinO 


n  B  sinu  m 

£  -S - ^  +  i  £  r.ta|r, 

1=1  ra  2:1  b=l  b  1 


(5.3) 


(5.-4) 


Noting: 


Rsine  =  r.sine,  +  risinai 

J  J  J  J 

Rcos0=ri  cos0i  +  r .  cosai 

U  J  J  J 

Tj2  =  R2  +  Rj2  -  2RRjCOs(0-0j ) 


One  may  write:  n 

(t>^^  =  RcosO  +  Z 


B  (Rcos"- ,R  cos  ) 
a  d  a  d  a 


a-l  R2+(d«a)2-2RdVOS<1-d'a) 

m  ,  r  Rcos’1-  R,  cos  0, 

-1  v  d  v  b 


1  -1 
-J-  Z  F,  cos 

2"  b=i  b 


_  vd  v "  b _ 

R2+(  R,  )  2-2R  R  cos(0-  0.  ) 
v  b  v  b  v  b 


(5.5) 


(i) 

ip  =  Rsin:*  +  Z 


B  (  R  sin,''  -Rsin  ) 
a  d  a  da 


a=l  R  +(‘ R  )  -2R  ,R  cos(  -  ) 

da  da  da 

1  m  ?  a 

+  77  £  F  :n  R-+(  R  “)-2R  R,  cos( —  ,)|  (5.6) 


where  the  subscripts  d  and  v  denote  quantities  due  to  doublets  and 
vortices,  respectively. 

Now  3a>  dRg,  dea,  rb,  vRb,  and  v0b  must  be  determined. 
For  n  doublets  and  m  vortices,  there  are  3(n+m)  unknowns.  Thus, 
specifying  3(n+m)  coorindate  points  (R,0)  of  the  profile  =  0  , 

results  in  a  system  of  (3n+m)  nonlinear  equations  in  3(n+m) 
unknowns.  Such  a  system  may  be  solved  by  a  nonlinear  counterpart  of 
the  Gauss-Seidel  iteration  method. 

If  3(n+m)  coordinate  points  of  <|/  ^  =  0  are  known,  the 
one  may  set  up  the  following  system: 

0  =  f1 (X1 ,X2,  .  .  .  ,  ^3(n+m)) 

0  =  f2(,-x1  ’x2>  *  *  *  *  X3(n+m)  ^ 


0  f3(n+m) (X1 ,X2’  ’  • 

Letting: 


■1 

Xn 

be 

n+1 

"*  Xn+m 

be 

n+m+1 

X2n+m 

be 

2n+m+1 

X2(n+m) 

be 

2(n+m)+1 

"*’X2(n+m) 

be 

3n+2m+ 1 

X3(n+m) 

be 

X3(n+m) 5 


d  1 


^dRn 


vR1  %Rm 


d  1 


~*den 


v®1  >  v9n 


and  1 et t i ng  superscr i p ts  k  and  k  +  1 


denote  successive  iteration 


steps,  superscript  k*  denote  the  value  after  the  kth  step  unless 
the  value  of  step  k  +  1  is  available,  and  subscript  l  denote  the 

given  coordinate  point  of  1/  i  ^  =  0  .  The  following  iteration 
equations  result.  See  following  pages. 

There  are  two  healthy  advantages  to  this  Gauss-Seidel 
iteration  method.  First,  truncation  and  round-off  errors  do  not 
accumulate  as  they  do  in  elimination  methods.  Each  new  value  is,  in  a 
sense,  a  new  initial  guess.  Second,  a  relatively  small  amount  of 
memory  is  required.  With  N  total  unknown  and  N  specified 
coordinate  points,  only  2N  memory  locations  are  needed.  Contrast 
this  to  some  elimination  methods  which  require  N2  memory  locations. 

However,  due  to  the  strongly  non-linear  nature  of  the 
equations  to  be  solved,  one  must  carefully  formulate  an  initial  guess. 
The  following  procedure  is  suggested: 

1.  Build  a  symmetrical  profile  with  a  thickness  distribution 
roughly  equivalent  to  the  desired  non-symmetrical  profile. 
Thus,  fj  =  vRj  =  vBj  =  d9j  =  0  . 

2.  Solve  the  first  and  third  equations  of  the  above  iteration 
equations  subject  to  step  1.  Make  initial  guesses: 

Bj°  =  (Thickness  at  j)/n 
dRj°  =  (j/n)  (chord  length) 

3.  Change  the  symmetrical  profile  to  a  slightly  non-symmetrical 
one. 

iJ.  Use  the  Bj  and  dRj  found  by  step  2  as  an  initial  guess 
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^  B  (,R  sin -R.sir. 
a  d  a  d  a 

r2+(  Rk)2-2R.  Rkco.s (•■■,-  k) 

-1  da  vd  a  *.  d  a 


^  ^  k  *  .  9  b  '>  b 

>  ",  ■  n  I  R?+(  R,  )  *"- 2 R  .  R,  cos) 


along  with 

dej  ■  »6j  •  0 

vRj  =  (j/m)  (chord  length) 

Tj°  =  .01  sin(j-rr/m) 

5.  Solve  the  above  iteration  equations  subject  to  step  4. 

6.  Repeat  steps  3  through  5  until  the  desired  profile  is 
obtained. 

During  the  initial  design  stages  one  is  interested  in 
finding  a  rough  profile  with  desirable  flow  characteristics.  It  is 
economical,  therefore,  to  specify  fewer  coordinate  points  in  this 
stage  than  in  later  stages  where  the  precise  profile  is  required.  The 
above  procedure  need  only  be  employed  in  the  initial  stage,  however. 
As  one  specifies  more  coordinate  points,  one  will  have  to  specify 
additional  initial  guesses  for  additional  doublets  and  vortices. 
These  additional  initial  guesses  may  be  intelligently  chosen  by 
interpolation  of  the  final  values  yielded  by  the  above  procedure. 

Leading  and  trailing  edge  closure  will  always  be  a  problem 
with  a  finite  number  of  doublets  and  vortices.  Luckily,  the  non¬ 
closures  may  be  reduced  to  negligible  magnitudes  by  specifying  closely 
spaced  coordinate  points  at  the  nose  and  tail  of  the  profile.  The 
non-closure  is  thereby  forced  to  fall  between  a  particular  pair  of 
arbitrarily  close  points. 

The  angle  of  attack  is  easily  found  as  the  angle  between  the 
chord  line  and  the  normal  to  the  profile  at  the  front  stagnation 


The  velocity  field  data  are  now  at  hand.  Specification  of 
any  two  points  in  the  field  allows  the  calculation  of  A<J>^  and 
Aij/^  by  equations  (5.5)  and  (5.6).  If  the  two  points  are 
sufficiently  close  to  each  other,  then 

A<f/  *^=d<j>^ * udx  +  vdy  (5.13) 

Atj/^=dif^^=  -vdx+udy  (5.1*0 

Solution  of  equations  (5.13)  and  (5.1*1)  yields  u  and  v  .  By 

q  =v/u^+v2  (5.15) 

9  =  tan  ^v/u)  (5.16) 

The  required  data  are  known  throughout  the  field. 


Chapter  6 


DISCUSSION 

The  hodograph  method  presented  in  Chapter  4  extends 
Lighthill's  transformation  method,  valid  for  subsonic  flows  only,  to 
transonic  flows.  The  ability  to  treat  supersonic  regions  springs  from 
theorem  4. 

An  underlying  assumption  in  theorem  4  is  that  t  is  a 
+  ^  t 

function  of  e~  .  This  assumption  predicts  that  x  and  t  vary 
together,  i.e.  they  both  increase  or  they  both  decrease.  Lighthill, 
however,  assumed  x  to  be  a  function  of  e±2t  in  his  development  of 
the  asymptotic  formulae  of  the  Chaplygin  functions  for  supersonic 
flow.  This  assumption  is  physically  insensible  as  it  predicts  x  is 
constant  in  supersonic  flow.  Hence,  the  ratio  (q/qm)  is  constant 
irrespective  of  t  . 

Consulting  Lighthill's^  work  shows  the  form  of  theorem  4  is 
similar  to  Lighthill's  analogous  result.  Both  predict  poles  at 
negative  integral  n  as  well  as  the  oscillatory  behavior  of  ^n(x) 
along  the  positive  real  axis  when  x  >  x.  .  Both  results  predict  that 
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for  large,  positive,  real  n  the  zeroes  of  ( t )  are  given 
approximately  by  the  zeroes  of  the  trigonometrical  functions.  Theorem 
4,  however,  is  far  less  complicated  than  Lighthill's  analogous 
formulae.  From  this  simplicity  rises  the  transformation  technique 
developed  in  this  study  which  is  markedly  similar  to  Lighthill's 
method  but  is  able  to  treat  transonic  flows. 

All  the  advantages  of  Lighthill's  technique  are  possessed  by 
this  newly  developed  method.  In  particular,  any  flow  which  can  be 
represented  by  potential  and  stream  functions  may  be  transformed 
utilizing  purely  numerical  data  concerning  the  flow  velocity  magnitude 
and  direction.  No  i.onger  must  one  be  restricted  to  those  flows  which 
can  be  represented  as  Laurent  series  expansions  or  Mel  1  in-Barnes 
integrals  in  the  hodograph  plane.  The  governing  incompressible 
complex  potential  about  a  general  lifting  profile  is  generally  far  too 
complicated  to  be  represented  thusly  and  cannot  be  employed  as  an 
incompressible  boundary  condition.  The  new  method,  however,  is  able 
to  transform  such  a  boundary  condition. 

This  method  eliminates  the  closure  problem  experienced  by 
those  who  use  such  representations  of  the  model  flow  in  the  hodograph 
plane.  Previous  researchers  using  hodograph  methods  to  solve  the 
compressible  transonic  flow  about  a  profile  have  overlooked  another 
work  due  to  Lighthill11.  Lighthill  addressed  the  problem  of  inverse 
design  and  formulated  three  constraints  which  must  be  satisfied  if  the 
inverse  problem  is  to  be  well-posed.  They  are: 


where  q  (m)  is  the  prescribed  speed  distribution  as  a  function  of  the 
polar  angle  of  the  circle  obtained  from  a  conformal  mapping  of  the 
airfoil  in  question.  Qualitatively,  the  first  constraint  states  that 
the  pressure  distribution  over  the  airfoil  and  the  free  stream  speed 
cannot  be  specified  independently  of  each  other.  The  other  two  state 
that  the  specified  pressure  distribution  and  the  angle  of  attack  may 
not  be  independently  specified  and  that  the  trailing  edge  close. 
Trailing  edge  closure  problems  are,  therefore,  symptomatic  of  ill- 
posed  inverse  problems.  Those  researchers  using  Laurent  series 
expansions  and  Mel  1  in-Barnes  integrals  to  represent  their  model  flow 
in  the  hodograph  plane  are  specifying  a  pressure  distribution  around 
some  profile.  They  must  arbitrarily  assume  a  free  stream  speed  to 
begin  their  solution  process.  In  so  doing,  they  violate  Lighthill's 
first  constraint  and  closure  problems  result.  Given  a  closed  profile 
in  incompressible  flow,  the  new  transformation  method  will  transform 
it  into  another  closed  profile. 

The  construction  of  a  model  flow  about  a  closed  lifting 
profile  is  only  theoretically  possible  by  the  method  of  Chapter  5.  A 
closed  lifting  profile  can  result  only  when  the  number  of  vortices  is 
infinite.  Thus,  in  practice  a  small  closure  problem  will  exist  over 
some  small  interval  at  the  trailing  edge  and/or  the  leading  edge.  By 
specifying  the  coordinate  points  at  the  nose  and  trailing  edge  of  the 
given  profile  arbitrarily  close  together,  the  closure  problem  is 
forced  to  fall  over  some  arbitrarily  small  interval.  Thus,  the 


closure  problem  can  be  reduced  to  effects  which  may  be  ignored  for 
engineering  purposes. 

Another  possible  drawback  to  the  construction  of  the  model 
flow  by  the  method  presented  in  Chapter  5  lies  with  the  initial  guess 
of  the  model  flow.  The  stability  of  the  numerical  method  is  strongly 
dependent  upon  this  initial  guess.  However,  if  the  procedure  outlined 
in  Chapter  5  is  followed,  a  reasonably  accurate  initial  guess  may  be 
constructed  which  should  eliminate  the  stability  problem. 

The  limiting  maximum  local  Mach  number  of  2.2735  may  be  a 
theoretical  drawback,  but  not  a  practical  one.  The  highest  local  Mach 
number  existing  in  shock-free  flow  about  any  airfoil  this  author  has 

seen  is  1.42.  This  Mach  number  existed  over  two  airfoils  designed  by 
12 


Boerstoel 


Chapter  7 

CONCLUSION 

M.  James  Lighthill  in  19^7  developed  a  method  of  solution  of 
the  hodograph  equation  governing  the  subsonic  flow  around  a  body. 
Employing  the  incompressible  flow  around  a  body  as  a  boundary 
condition,  this  method  transforms  the  boundary  condition  into  a 
solution  of  the  hodograph  equation.  Preventing  such  a  transformation 
in  supersonic  flow  are  the  complicated  and  physically  insensible 
supersonic  asymptotic  formulae  of  the  Chaplygin  functions  which 
Lighthill  developed. 

This  study  developed  a  phys i ca 1 1 yr easonab 1 e  asymptotic 
formula  of  the  Chaplygin  functions  valid  for  supersonic  flow  up  to 
Mach  2.2735  in  air.  The  new  formula  permits  the  development  of  a 
method  which  transforms  an  incompressible  boundary  condition  into  a 
solution  of  the  hodograph  equation  governing  both  subsonic  and 
supersonic  flow  about  a  body.  Any  incompressible  flow  around  a 
profile  that  is  governed  by  an  incompressible  complex  potential  may  be 
transformed  into  a  solution  representing  the  compressible,  transonic. 


two-dimensional,  potential  flow  about  a  similar  profile.  This 
transformation  is  developed  from,  and  is  markedly  similar  to, 
Lighthill's  transformat  ion. 

Only  numerical  data  concerning  the  flow  velocity  magnitude 
and  direction  around  a  given  profile  are  required  for  the 
transformation.  Numerical  methods  required  by  this  method  are  simple 
and  well  behaved. 
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Appendix  A 

DERIVATION  OF  THE  GOVERNING  EQUATIONS 


Any  study  in  fluid  dynamics  must  begin  with  the  governing 
equation(s)  of  motion.  The  governing  equation  is  often  derived  under 
certain  assumptions  which  limit  its  use.  These  assumptions  must  be 
known  so  that  the  governing  equation(s)  may  not  be  incorrectly  applied 
to  certain  flows. 

The  following  derivations  parallel  those  of  Scnapiro1. 

Beginning  with  the  assumption  of  irrotational  motion  it  is 
easily  shown  that  the  integral 


B 


J  q  cosadil 
A 


(  A.  1  ) 


is  independent  of  path.  Thus,  q  cosadS.  is  an  exact  differential  and 


one  may  set: 

d<j>  =  q  cosadt 

where  <p  is  the  velocity  potential, 
flow  are  often  used  interchangeably. 
ipx  =  d<p/dx  =  u 


(A. 2) 

Potential  flow  and  irrotational 
Thus,  for  two-dimensional  flow 

(A. 3) 


5h 


3$/3y  =  v 


(A.  t) 


where  u  and  v  are  the  velocity  components  in  the  x- 


and  y- 


directions,  respectively,  of  a  Cartesian  coordinate  system.  The  above 
is  summarized  as: 

q  =  V4>  (A. 5) 

For  a  steady  two-dimensional  flow,  the  continuity  equation 

becomes 

3(pu)/3x  +  3(pv)/3y  =  0  .  (A. 6) 

Carrying  out  the  indicated  differentions  yields 


p^xx+<*Vy)  +  4>x3p/3x  +  <J>y3p/3y  =  0  . 


( A. 7) 

Next,  an  equation  is  required  relating  p  with  the  potential 
function.  Summing  Euler's  equations  of  motion  along  a  streamline  in 
steady  flow  yields: 


dp  =  -pd(q^/2; 


Expanding: 


2  2 

dp  =  _  ,d(^4^)  =  -od(J5>)  . 


2  '  2 
For  isentropic  flow,  the  sound  velocity  is 

2 


dp/dp 


Substituting  into  the  above 


3x 


A  =  ~  —  (dp)  -  -  ~  ( .  '  +  :  :  ) 

Jx  2  2  ix  2  x  xx  v  vx 

c  c  c  ' 


Likewise: 

3p/3y  =  -p/ c^( <J>x4>xy  +  4>yyyy )  . 
Substitution  of  the  above  into  equation  (A. 7)  yields: 


(A. 3) 


(A. 9) 


(A.  10) 


(A. 11 ) 


(  A. 1 7) 


The  local  sonic  speed  is  derived  from  the  energy  equation  for  an 
isentropic  flow  of  an  ideal  gas.  It  is: 

2  2  Y-l  2  2  y-l  ,  2  2 

c  =c  -  ~ g—  q  =c  -  -or- 

o  2  o  2  xy  ( A  #  1 4 ) 

Note  that  the  isentrope  for  an  ideal  gas  is  given  by 

p/pY  =  constant  (A. 16) 

from  which  one  can  see  that  specifying  a  straight  line  isentrope  in 
the  p  vs.  1/p  plane  amounts  to  specifying  Y  =  -1  . 

Equation  (3)  is  the  governing  equation  of  the  potential 
function.  Obviously  non-linear,  its  solution  is  at  best  difficult 
unless  simplifications  are  made  or  it  is  somehow  transformed  to  a 
linear  equation. 

The  stream  function  is  now  defined.  The  stream  function,  \|>  , 
exists  only  for  the  two-dimensional  steady  flows  and  is  defined  by: 

<PX  =  U  ~  Pc/p^y  1  4>y  =  V  =  -pQ/pipx  .  (A.  16) 

Substituting  the  above  into  equations  (A. 13)  and  (A. 14)  yields  the 
governing  partial  differential  equation  of  the  stream  function: 


For  an  ideal  ga3 


CpT  =  ^Cp/YR)  YRT  =  ( 1 /Y~1 ) . 


(A. 20) 


Combining  these  two  equations  and  evaluating  the  constant  at  three 
reference  conditions  (zero  speed,  zero  temperature,  and  sonic  speed) 
yields: 

q2  +  [2/Y-l]  c2  =[2/Y-l]  eQ2  (A.  21) 

=  Qm2  (A. 22) 

=  (Y+1/Y-1)  c2  =  (Y+1/Y-1)  q2  .  (A. 23) 

Equations  (A. 21)  and  (A. 22)  may  be  recast  as 

2  p 

J_  P  +  =  _I_  11  (A. 24) 

Y-l  P  2  y-1  o 

(2Y/Y-1)  pQ/p0  =  qm2  (A. 25) 

p 

by  noting  c^  =  Yp/p  for  isentropic  flow.  Solving  equation  (A. 24) 
for  q  and  employing  the  isentrope  relation  yields  the  De  Saint 

Venant  and  Wantzel  formula: 

2  2 , .  ,  :  .  v '-1 , 

q  =  qmU  ”  26) 

o 

The  final  required  relation  is  found  by  substituting  the  above  into 
equation  (A. 21)  yielding: 


Appendix  3 


The  Molenbroek-Chaplygin  Transformation 

In  equations  ( A.  1 3 )  and  (A.17)  the  independent  variables  are 
the  Cartesian  coordinates  (x,y).  By  considering  these  coordinates  as 
functions  of  the  polar  velocity  coordinates  (q,0)  the  governing 
equations  may  be  linearized.  This  is  the  Molenbroek-Chaplygin 
transformation. 

The  following  derivations  follow  those  of  Boerstoei2  and 

O 

Ferrari  and  Tricomi  . 

It  follows  from  the  definitions  of  the  stream  and  potential 
functions  that  their  total  derivations  may  be  related  thusly: 


d<J>+i(p0/p)d^=udx+vdy+i(udy-vdx) 

Thus: 

=(u~i v) (dx+ idy ) =qe  l9dz 

(3.  1  ) 

dz=ei  9/q[dq+i(p0p)dt{/] 

(3.2) 

3z/9q=ei0/q[<i>q  +  i(po/pHcl] 

(3.3) 

3z/30=ei  9/q[<J>g  +  i  (pQ/p)  <|>Q] 

(3.«0 

When  x 

and  y  are  considered  as  continuously 

differentiable 

functions  of  q  and  0  ,  their  mixed  partial  derivatives  exist  and 
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the  order  or  differentiation  is  immaterial.  Hence 


{e19/q[<{>q  +  i(p0/p)^Q]}9 


=  {e10[<f>0  +  i(po/ p)tp0]}q 


yielding: 


ie19/qUq  +  i(p0/p)ij/q] 

=  -el9/q2  46+iel3(p0/pq)q^9  . 
Equating  the  real  and  imaginary  parts  yields: 
<i>e=q(Po/p)*q 


fq(p0/pq)q'J'e 


( B.  6) 


(B.7) 

(B.3) 


We  define 


t=(q/qm)" 


( B.  9) 


Subtracting  equation  (A. 20)  from  equation  (A. 26)  yields: 


PQ/p=( 1-t) 


-1  /Y-1 


( B. 1 0) 


Equations  (B.7)  and  (B.8)  can  now  be  expressed  in  terms  of  x  . 
Noting 


dx/dq=2q/ q 


(B.11) 


we  obtain: 


<J>0  =  2(q/qm)2(po/P)^T 
4>0=2t(1-t)"1/Y"1^t 


Q=2t(1-t) 


•1  /Y-1 


(B.12) 


(B.13) 
(B.1U) 
( B. 1 5) 


Likewise: 


bT=q(po/pq)T^0 

3ivoq  q  di  .•  di  q 

q  3/3x(pn/pq)=( Y-1 )"1(1-t)_Y/Y"1-1/2t(1-t)'1 /Y_1 


( B. 1 6) 


(B.1 7) 


(B. 1 3) 


-1/2t(1-t)  "  '  ' L 1 -T- ( 2t/Y-1 ) ] 
-1/2t(1-t)Y/Y_1[1-t(Y+1/Y-1)]  . 


( B.  1  9 


( B. 20) 


From  equation  (A. 20) 


l2/2  +  c2/Y-1  =  q  2/2 


we  see  when  q  =  c 


Hence: 


T  _  =  ( c/q_ ) 2  =  Y-1/Y+1  . 


*0  =  P*T 

P=  -2t(1-t)Y/Y  Vd-T/Tg) 


( B. 21  ) 


( B. 22) 
(B.23) 


Cross-differentiating  equations  ( B.  1  -4 )  and  B.22)  to  eliminate  the 


potential  functions  yields: 

PQ*tt+PQtV’J’90=s° 


(B.2H) 


This  is  the  linear,  governing  mixed  partial  differential  equation  of 
motion  in  the  hodograph  plane  in  terms  of  the  stream  function.  The 


discriminant  of  the  above  is: 


B2-iJAC=0-4(PQ)(-1  )=JJPQ 


(B. 25) 


The  coefficient  PQ  is  negative  when  t<ts  .  Thus,  equation  (B.214) 
is  elliptic  for  subsonic  flow.  Likewise,  PQ  is  positive  when  t>ts 
and  equation  (B.24)  is  hyperbolic  for  supersonic  flow.  The  mixed 
behavior  is  evident. 

The  variables  t  and  0  are  separated  by 
H>(t  9)=^n(T)e±ine 


yielding: 


PQ*nTT+PVnT+n  V° 


( B. 26) 


Substituting  ^n(i)  =  Tn  "Fn(i)  into  equation  (B..?6)  yields: 


'  „  ‘  »  *  *  .  »  -  . 


PQF»+[PQ(n/t2)+PQT]F^ 

+[PQ(n/2) (n/2~1 )x'2+PQ(n/2x)+n2]Fn=0  (B.27) 

where  the  prime  denotes  differentiation  with  respect  to  x  . 
Substituting: 

PQ=  -4t2(1-t)/(1-t/ts)  ; 

PQt=  -i»T[1+(2-Y)/(Y-1)x]/(1-T/Ts)  (B.  28) 

into  equation  (B.27),  collecting  and  cancelling  terms,  yields: 
xO-T)F"(x)+[n+1-(n-(2-Y)/(Y-1)T]F^(T) 

+n(n+1 )/2( Y-1 )  Fn(x)  =  0  (B.29) 

Comparison  with  the  hypergeometric  equation  of  which  the 
Gaussian  hypergeometric  functions  are  solutions  shows  F  (x)  is  a 
Gaussian  hypergeometric  function.  In  this  case: 

anbn  =  -n(n+1)/2(Y-1)  (B.30) 

an+bn  =  n-(1/Y-1 )  (B.31 ) 

cn  =  n+1  .  ( B. 32) 

Solving  for  an,bn  yields: 

an,bn=1/2[n-1/Y-1±((Y+1)n2/Y-1  -  (1 /Y-1 )2) 1 /2]  ( B. 33) 

where,  by  convention,  aR  >  bn  .  The  Gaussian  hypergeometric 
functions,  found  by  the  Method  of  Frobenius  about  x  =  0  (other 
regular  singular  points  are  1,  °°  ),  are: 

F(a,b;c;x)  =  1 *( ab/c) x+a(a+1 )b (b+1 ) /2 ! c( c+1 ) x2 

+(a(a+1 ) (a+2) (b) (b+l ) (b+2) /3 ! c(c+1 ) (e+2) 

(B.3'0 

The  Chaplygin  functions,  ^n(x)  ,  are  defined  as: 

^n(x)=xn/2F(an,bn;n+1 ;x) 


(B.35) 


Correspondence  between  the  hodograph  and  physical  planes  is 
one-to-one  when  the  Jacobian  of  the  transformation  is  not  zero  or 
infinite. 

J  =  9(x,y)/3(q, 0) 

=  (3)/3q)(3y/30)-(3x/3e)(3y/3q)  (B.36) 

From  equations  ( B.  3 )  and  (B.4),  and  noting  Z  =  x  +  iy  ,  the  Jacobian  is 
seen  to  be 

J  =  e2l9/q2  (  pQ/ p )  C H  •  ( B.  37) 

Substituting  equations  (B.7)  and  (B.8)  into  the  above  yields: 

J  =  e2l9/q(po/p)[(po/pq)qij;02-(po/p)\J;q2]  ( B.  38) 

Noting: 

3/3q(pQ/pq;  =( pQ/p ) (-q~2 ) +( pQ/q)d/dq(  1 /p) 

=  -p0/pq2-^p0/pc2 

=  pQ/p(1  /c2-1  /q2)  (B.39) 

equation  ( B . 3 8 )  becomes: 

J  =  e2i9/q(p0/p)2[(1/c2-1/q2)4;2-^2]  (B.40) 

J  =  -e2i0/q3(po/p2)[q2ij>q2-(M2-1  H02]  (B.JJ1) 

For  a  one-to-one  correspondence  between  the  hodograph  and  physical 
planes,  we  must  have: 

q*  o;  q2^2-(M2-1  HQ2*o,“  (B.^42) 

When  the  Jacobian  is  infinite,  branch  points  or  branch  lines 
occur.  Branch-points  exist  in  subsonic  flow  and  branch  lines  exist  in 
supersonic  flow.  An  obvious  branch  point  is  the  stagnation  point 
where  q  =  o  .  Ferrari  and  Tricomi'  ,  show  that  no  singularity  arises 
in  the  fluid  flow  at  branch  line  images  for  supersonic  flow. 


6 


Limit  lines  occur  when  the  Jacobian  is  zero.  In  subsonic 
flows,  the  Jacobian  could  be  zero  at  the  stagnation  point  or  at  a 
point  of  infinite  curvature.  Supersonic  flows,  however,  can  have 
lines  where  the  Jacobian  is  zero.  The  equations  defining  these  lines 
are: 

q<i»q  <J>0  =  o 

q^q  +/m2-1  <1j0  =  o 

A  theorem  due  to  von  Karman  states  that  in  the 
plane  the  limit  lines  are  the  loci  of  points  of  tangency  of 
lines  with  the  characteristics  of  the  governing  equation, 
of  a  characteristic  in  this  case  is: 

de/dx  =  ~B±  \f  B2-1!  AC  /2 A  =  +  1 /^/“PQ 
See  Figure  2. 

A  qualitative  discussion  will  help  the  visualization  of 
limit  lines.  Because  the  Jacobian  is  zero,  a  finite  length  in  the 
hodograph  plane  is  mapped  to  zero  length  in  the  physical  plane.  Thus, 
the  physical  plane  could  be  considered  to  be  folded  at  the  image  of  a 
limit  line.  Although  the  velocity  will  be  a  smooth  function  about 
the  limit  line  in  the  hodograph,  it  will  be  discontinuous  in  the 
physical  plane  at  the  limit  line  image.  The  presence  of  a  limit  line 


(8. 43) 
(B.44) 
hodograph 
the  stream 
The  slope 

(B.H5) 


predicts  a  shock  in  the  physical  flow. 


Appendix  C 


THE  HODOGRAPH  PLANE  FOR 
A  TYPICAL  LIFTING  AIRFOIL 

The  nature  of  the  hodograph  plane  depicting  the  flow  about  a 
typical  lifting  airfoil  will  be  illuminated  by  three  examples. 

Non-Lifting  Circular  Cylinder  in  Incompressible  Flow 

Consider  first  the  complex  potential,  4>^ i  ^  +  i^*)  , 


of  an  incompressible  uniform  flow  about  a  circle.  With  the  circle 
radius  and  free  stream  velocity  normalized,  this  potential  is 


*(i) 

=  Z  +  Z„  1 
c  c 

(C.1  ) 

and 

d$(  1 } 

/dzc  =  u-iv  =  qei9  =  1-ZC2 

(C.2) 

Defining 

i  ft 

qe  and  substituting  into 

the 

above  yields: 

Zc  " 

±[o-ccr1/2 

( C .  3 ) 

=  ±[d-tc)1/2  -  o-ccr1/2] 

• 

(  C  .  U  ) 

Thus, 

tc  must  be  defined  on  a 

two- 

sheeted  Riemann 

surface 

for  an  unequivocal  one-to-one  correspondence  to  exist.  French  points 
exist  when 

dCc/dZc  =  0»°°  •  <r  ■■  ' 

Performing  the  differentiation  yields 
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d?c/dZc  =  2/Zc3  (C.6) 

Possible  branch  points  exist  when  ZQ  =  o  and  when  ZQ  =  * 

This  point,  not  in  the  external  flow,  is  therefore  of  no  interest. 
When  ZQ  =  «  ,  5=1  and  one  is  at  the  force  stream  conditions. 
The  branch  cut  extends  from  =  1  to  infinity  along  the  real  axis. 
Note  also  that  q  =  1  is  a  singular  point  of  equations  (C.3)  and 
(C.4). 

Consider  the  stream  function  about  the  circle.  It  will 
obviously  be  antisymmetric  with  respect  to  the  axis  of  flow.  Let's 
look,  therefore,  at  the  flow  above  the  axis  of  symmetry.  The  front 
stagnation  point  maps  to  the  point  (o  ,  tt/2)  on  the  hodograph.  At  the 
point  of  greatest  cross-section  presented  to  the  flow,  the  flow  angle 
is  zero  while  the  flow  magnitude  is  a  maximum.  Beyond  this  point,  the 
flow  angle  becomes  negative  and  the  magnitude  decreases  until  at  the 
rear  stagnation  point  the  velocity  is  zero  again  while  the  flow  angle 
is  -tt/2  .  See  Figure  3. 

Lifting  Circular  Cy 1 inder  i n  Incompressible  Flow 

Consider  next  the  same  situation  with  circulation.  The 
complex  potential  for  the  incompressi ble  flow  is  now: 

$(i)  =  Zc  ♦  Zc_1  *  i  f/pTT  £nZ(,  (C .  7) 

and 

d*(i)/dZc  -  u-iv  =  C0  =  1-ZC‘°  T  ijv.li,  zc_1  (C.3) 

We  see  the  hodograph  plane  is  again  two  sheeted.  Finding 
the  possible  branch  points: 
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(C.9) 


d^c/dZc  = 


ir/2ir  z„ 


Thus: 


dVdZc 

=  0 

at 

Zc  =  » 

( C .  1  0 ; 

dVdZc 

=  0 

at 

xc  =  -^i/r 

(C.  V  1 

d?c/dZc 

=S  CO 

at 

Zc  =  0 

(C.12) 

The  corresponding  points  in  the  hodograph  plane  are,  respectively: 


5C  =  1 

Cc  =  1  -  (T/iJir)2 

Cc=  " 


(  C  .  ;  ■}) 
(C.1  H) 
(CJ5X 


The  point  ZQ  =  o  is  again  of  no  practical  interest.  From 
equation  (C.8),  one  sees  that  5  =  1  can  be  reached  from  two  values 
of  Zp  ,  namely: 


Zc  -  ( C .  1  6 ) 

Zc  =  -2tt  i/r  (C.  1  7) 

Therefore,  the  point  5  =  1  cannot  be  the  branch  point.  Indeed,  the 

value  of  t  corresponding  to  the  free  stream  conditions  (Z„  =  »)  is 
a  singular  point  again  and  must  not  be  defined  on  the  sheet  containing 
the  regular  point  Z  =  -  2ni/r  .  One  is  left  with 

C*  =  1  ~  (  T/'Mtt ) 2  (C.1S) 

as  the  last  possibility.  Substituting  this  into  equation  (C.S)  shows 
that  5*  can  be  reached  from  one  and  only  one  point  in  the  oh  vs i cal 
plane.  Thus,  is  the  branch  point  and  the  brar.cn  cut  rune  along 

the  real  axis  out  to  infinity. 

Qualitatively,  the  physical  flow  image  in  t  j-.e  b.c.l.n-r  irh 
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plane  is  similar  to  the  non-circulatory  case.  Here,  however,  t Y 
stream  function  is  not  antisymmetric  with  respect  to  the  flow  axis  ~ 
the  maximum  velocity  magnitude  over  the  top  half  of  the  cylinder  i 
greater  than  that  obtained  over  the  lower  half.  See  Figure  t. 


Lifting  Airfoil  in  Compressible  Flow 

Consider  a  lifting  airfoil.  The  incompressible  oomple 
potential  may  be  conformally  mapped  from  an  incompressible  flo 
representing  that  about  a  circle.  This  general  mapping  i 


( C-lauert1  3 ) : 


Z  =  7. 


(C.1? 


where  the  coefficients,  ai  ,  are  generally  complex.  Now: 

r  -  dl>(i)/dz  =  (d0(  U/dZc)(dZc/dZ)  =  cc(dZ/dZc)_1 
dZ/dZc  =  1  -  aiZc'2-2apZc"3  ...  . 


(C.  21 


Assume  ZQ  larg  enough  to  neglect  orders  greater  than 


(dZ/dZ, 


=  1+aiZc-2+ai2Zc-^ 


(C.<?2 


Equation  (C.20)  becomes: 


e  =  (1-2C  2  +  ir/2TT  Zc  1)(i^a1zc  2+3l 2; 
C  =  1  +ir/2iT  -(1-ai)Zc~2+0(Zc-3) 


( C . 


Comparing  equations  (C.8)  and  (0.2^)  shows  the  lifting  airfoil  i: 
similar  to  the  lifting  circular  cylinder  except  that  the  branch  poin' 
has  moved  off  the  real  axis  and  the  branch  cut  is  now  a  radial  ra 
beginning  at  the  branch  point 


r*  =  1  -f/l  6tt‘-  ( 1  -a, ) 


/■•■nding  to  infinity.  The  point,  corr "r.pondi nr  to  the  free 


conditions,  5  =  1  ,  is  still  a  singular  point. 

The  'nodograph  plane  still  has  the  same  basic  structure  in 
compressible  flow  as  it  did  in  incompressible  flow.  For  low  speeds 
one  may  use  small  perturbation  theory  to  linearize  the  governing 
equation,  (A. 13).  The  result  is: 

32<{>/dx2  +  32/<j)/3y2  =  0 

6  co  = 

Through  Gothert's  Rule  and  the  Prandtl-Glauert  Similarity  Rules,  one 
can  predict  that  the  compressible  complex  potential  is  the  same  as  the 
i ncompr ess i b 1 e  case  for  affinely  related  airfoils.  However,  the 
previous  result  is  for  a  perfectly  general  airfoil  in  incompressible 
flow.  Thus,  the  same  result  will  hold  for  any  general  airfoil 
affinely  related  to  this  general  airfoil.  In  particular,  the  point  t, 
=  1  is  still  a  singular  point  and  the  hodograph  is  still  two-sheeted. 


See  Figure 
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Appendix  D 


INTRODUCTORY  RELATIONS  REQUIRED  IN 
ASYMPTOTIC  FORMULAE  DERIVATIONS 


Recall  the  characteristic  equations  of  equation  (1.6): 
de/dx  =  ±  (PQ)-1 72  (D. 1 ) 


Consider  first  the  case  where  de/dx  is  real,  i.e.  supersonic  flow. 
Then, 


de/dx  =  ±  1/2x((x/xs-1)/(1-x)) 


1  /2 


Hence: 


d0/dx  =  ±  1 /2x([(Y+1 )x-(Y-1 ) ]/(1-x)(Y-1 J)1 72  . 

/t-  j*  T 

de  =  ±/((Y+1  )x-(Y-1)/(1-x)(Y-1))172  dx/2x 

-'Y-l  /Y+1 


( D .  2) 
(D.3) 


(D.H) 


t  = 


tan 


-1  /(,+.!)  i-(-.-l) 


(l-0(,+D 


-tan 


-1  /(,+!) 


(!-■ )(.-!) 


(D.  5) 


See  Figure  6  for  a  plot  of  t  versus  x  .  The  governing  equation 
(1.6)  for  strictly  supersonic  flow  becomes 

PQC^Vr^x^Wx^'eo--0 


(D.  6) 


t‘T  =  (PC) 


-1 


^tt+(PQtxx+PQxtx)^r^eo=0 


(D.  7) 
(D.3) 


or 


^,tt  +  S^t^t  ^00~° 

3(  t )  =PQ[  +  (  PQ) _^72  (  P  Q+Q  )  ]  +  PQ  (  ±(  PQ)"1  72  ) 


(D.9) 


(D.  10 


S( t ) =±1 /?(  PQ)  ,/-(PQt-?tQ) 

S(t)=±1 /2(Qt/Q-Pt/P) 

=±2( Y  +  1 )/(Y-1 )2x2( 1 -t)_1 /2( ( Y  +  1  /Y-1  )t-1 )"3/2 


(D. 11  ) 


(D.1?) 


In  subsonic  flow,  ?  <  0  and  the  characterist i cs 


imaginary.  One  may  write: 


i  de/dx  =  ±i(PQ)  1/2  =  ±(-PQ)  1/2 


( D.  1  3 


Because  the  characteristics  are  imaginary  and  do  not  physically  exist, 
Lighthill-*  is  valid  in  redefining  the  character! Stic  equation  as 

de/dx  =  (-PQ)_1/2  .  (D.1U) 


Thus: 


f  d6  =  i  f ( ( ( Y+1  ) x-( Y+l) )/(1 -x) ( Y-1 ) ) 1 /2  dx/2x  . 

J  •'a 


(D. 15) 


Integrating, 


»-o+A± i  canlC 1  -H)_L  lilU 

y.-l  Y  (,  +  l)(l-)  Y  ('-1X1-:)  CD.  1 6) 

where  o  is  an  arbitrary  constant  yet  to  be  determined.  For  strictly 


subsonic  flow,  the  governing  equation  (1.6)  becomes: 

^ss  +  ^99  -  T(s^s  CD.  1  7) 

T(s)=-1/2(-PQ)1/2[Qt/Q-Pt/P] 

=2( Y  +  1 ) /( Y~1 ) 2  x2(1-x)"1/2(1-(Y+1/Y-1)x)_3/2  (D.18) 

Li ghthi 1 1 3  defines  o  to  b°  that  value  which  causes  x  to 
be  asymptotically  e2s  as  x— *>0  and  S — .  Symbolically,  this  is 
represented  as  x  ~  e  . 

This  exponential  form  is  particularly  well  suited  to  a 
series  solution  of  equation  (D.1 7 )  which  will  yield  the  important 
asymptotic  formula  for  the  Chaplygin  functions  for  subsonic  flow.  The 


point  t  =  o  ,  s  =  -  «  ,  at  which  t  ~  e  ,  represents  the 
incompressible  limit.  The  behavior  of  any  solution  to  compressible 
flow  -  which  will  include  ipn( t ) *s  -  must  reduce  to  the  solutions  of 
the  Laplace  equation  in  the  incompressi b 1 e  limit.  This  requires 
thorough  knowledge  of  the  behavior  of  t  as  t — *o  becomes  v„( t )  is 
a  function  t  .  This  matching  of  the  compressible  solutions  to  the 
incompressible  solutions  is  the  basis  of  Lighthill's  transform  ition 


method.  Thus, 

,-exp  i  *<« /$  -tamr1  . 


By  the  identity  tanh  1 ( x )  =  1/2  £njl+x/1-x|  ,  x2  <  1  ,  we  have: 


(D.19) 


.  2awl-x.  ,l+x\  V  Y-l 


(D. 20) 


Thus , 


>-!)-( .+l)r 

(>+l) (l-i) 


,  -  /(v-1  )-(>!) 

x  y  (,-Dd-o 


(D.21 ) 


t  ,  .  '  1+X;  /  .+1  :  1+X ' 


( D. 22) 


o=  ^njx  |-< 


As  t — *o  , 


St  /(  >-l  )-(  )  +  l  )  t-l  /'(  ,-l)-(,  +  D 

v-i  t""‘  y-feTui-.T  +t‘lnh  \/'T7-T)(i-V)- 


>-i )-(.+!) j 
(»-l)(l-v) 


(D. 23) 


(D.2U) 


The  second  and  third  terms  approach  -®  and  +»  ,  respectively. 


Using  the  binomial  expansion  and  neglecting  terms  of  second  order  and 


higher  in  t  ,  one  deduces: 


0.25) 


r t- 


r^nlrj 


+tanh 


Neglecting  the  x 


V  o-D (i-o 

term  yields: 


i  ;J2+  f(1- 

2 ' n  i  -1  j 

!  (,-D  ! 


1 /2  Hn|2( Y-1 ) |  .  (D. 26) 

Finally: 

o  =  1  /2Hn|2(  Y-1 )  |  -  Jy+T/Y--\  tanh"1  /yhVy+I"  .  (D.27) 

For  air,  Y  =  1.4  ,  o  =  -1.173  .  See  Figure  7  for  the  plot  of  s 
versus  x  .  Note  that  s — *-o  as  x — >xs  .  See  Figure  8  for  the 
plots  of  T(s)  and  +S(t)  versus  x  . 

Theorem  1  (Lighthill^)  If  o  <  x  -  1  ,  4in(x)  is  an 
analytic  function  of  n  except  at  n  =  -2, -3, -4,.  .  .  ,  where  it  has 
simple  poles,  its  residue  at  n  =  -m  being  -mCm^m(x)  ,  where  Cm  is 
positive  and  ~(2Ttm)_1  e~2om  as  ra— ►«  . 

For  any  m  , 

cm  =  r(am)r(1+m-bm)/r(an-m)r(1-bm)(m!)2  (D. 25) 

where  the  values  of  am  ,  bm  are  those  values  previously  defined  for 
the  hypergeometric  function,  equation  ( B . 33 )  • 

Lighthill-  also  deduced: 


4*0  =  1  (D.  29) 

0=  (Y-1){1-(1-x)Y/Y_1  }/Yy"x  (D.30) 

’I'i  =  x'1/2  +  (2Y-2)_1ii1  (x)  (D.31) 

Theorem  2  (Li  ghthill^)  If  <5  >  o  and  0^  <0  ,  then  Wn(s) — *1  ,  i.e. 
Vn(x)  ~  en°V(r)  ,  uniformly  for  s  -  o1  ,  and  for  n  in  the  whole 
complex  plane  with  circles  of  radius  <5  round  each  negative  integer 
excluded  as  [n| — . 


s 


Lighthill^  develops  theorem  7  by  assuming  a  form  of 


amenable  to  a  series 

qo  1 

ution  of  the  governing  equation  of 

suboor.i 

flow,  equation  (D.17). 

• 

First  integrate 

equation  CD. 1 8 ) 

f-P/Q  3 

1/22,n|P/Q| 

=  /  TCs,  ) ds. 

(D.  ?2 

1 

J  -*> 

r s 

1 /2£n | -P/Q | 

= 

/  T(s1)ds1 
**  —00 

f  D . 

(-P/Q)172  = 

EXP 

{  j  T( s1 )ds1 } 

J  —  00 

(  ^  Z  ;; 

Def ine, 

/s 

T(  s,  )ds1  }  ( D .  3 1 

-00 


where  the  additional  factor  of  1/2  is  needed  in  the  series  solution  of 
equation  (D.17).  Thus, 


V(o)  =  1  (D. 36) 

dV( t)/ds  =  1/2T(s)V(t)  .  (D.37) 

Lighthill^  assumes 

Yn(t)  =  ensV(t)Wn(s)  .  ( D . 3 S ) 

Note  that  ens  =  in/2  and  thus  V(i)Wn(s)  must  represent  F(an,  bn; 
n+1  ;  t)  .  Because  i|/n( t )  ,  e113  ,  V(t)  are  all  analytic  functions  of  t 
in  the  given  domain  (t*o  ,  ts)  ,  then  Wn(s)  must  also  be  an  analytic 

n  rv 

function  of  x  ,  and  hence  of  e °  .  Also,  when  x  =  o,  F  (  a  M  , 
bn;  n+1  ;  x)  =  1  .  Because  V(o)=  1  ,  so  too  must  Wn(-»)  =  1  .  Th:  s 
impl ies 


V3>  ■ 


"n.r'-' 


,T3 


r=o 


(D. 3?) 


Proof:  A s:'r-  it  true  for  m 


*r  <  it-r  )  Ac  ~  >4. 


-Ae  3  = 


'*r(n+r)  •  Ac 


A  -Jr  •  j  .  2rt  •  > 
3  - - L’  Cc  3 

+  D 


Now  4r(n  +  r) | >2 | n  |  5  ,  for  when  jr|>|n/2|  ic  i 
|r|<|n/2|  it  is  >  4jn/2|  .  Hence, 


••  !  •  ‘  >  -  ! 


,  -  2  r  •  t  ,  AC  ,  1.  Ac  1  - 

d  •  Ac  (  L+  - )  ~~  1  - 

n ,  r  ■  n .  2  n  n 


if  |n|  >A C/6  .  Hence  the  result  is  true  and  the  lemma  fellows, 


It  is  deduced  that,  for  s  -  o1  ,  under  the  restrict:  • 


1 emma  3 . 


Vs'"1  I  -  1  I  ^ n ,  r  I e"' 

r=  1 


(  D .  5 1  ) 


<  A/5 [ n |  ,  1/l-e"(o1  °2) 


and  thus  Wn(s) — *1  and  theorem  2  results.  One  would  do  well  to  note 
that  the  sum  after  K  terms,  SM  ,  of  the  series 


=  a  ( 1  -  r'  )  /  ( 1  -  r ) 


when  following  Lighthi ll’s  above  procedure.  Not' 


form  of  4r(n+r)  .  For  n  --  -r  the  above  is  i  n  v  a  1  :  .1  be  mice  |  d  |  is, 

indeterminate.  However,  ij;  ( t )  has  simple  poles  .at  t;  he'  negative  r'Vil 
integers  (except  n  =  -1)  and  Light  hi  l  1  is.  entir--'ly  valid  in  •'X'i  u.iir.g 


